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Systematics of Neotropical Spiny Mice, Genus 
Neacomys Thomas, 1900 (Rodentia: Cricetidae), 
from Southeastern Amazonia, with Descriptions of 
Three New Species 
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ABSTRACT 


Species of Neacomys are small cricetid rodents that occur in forested habitats of Central 
and South America, from eastern Panama to central Bolivia and central/western Brazil. In order 
to assess species diversity of this poorly known genus, we obtained cytochrome b gene sequences 
from the most comprehensive taxonomic and geographic sampling analyzed to date. We also 
conducted morphological analyses on a large series of specimens housed in 15 museums, 
including types of 10 out of 14 nominal taxa. Our analyses of the genetic data recovered 17 
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lineages clustered in four distinct clades. Among these lineages, 11 correspond to species cur- 
rently recognized as valid, and the remaining six are putative new species. In southeastern 
Amazonia—the geographical scope of this report—four undescribed species were discovered, 
three of which are named herein: Neacomys marajoara, sp. noy., from the Island of Marajo, Para 
state; Neacomys vossi, sp. nov., restricted to the Tapajds center of endemism (between the Tapa- 
jos and Xingu rivers); and Neacomys xingu, sp. nov., restricted to the Xingu center of endemism 
(between the Xingu and Araguaia/Tocantins rivers). The new species can be discriminated from 
other Neacomys species by the morphology of the nasal bones, zygomatic plate, interorbital 
region, subsquamosal fenestra, paraoccipital process, incisive foramina, auditory bullae, antero- 
cone and anteroloph of the first upper molar, carotid circulation pattern, and karyotype. Our 
results substantially improve our understanding of the genus Neacomys by providing morpho- 
logical, morphometric, and novel molecular insights about these poorly known rodents and 
demonstrate that the diversity of small Amazonian mammals is still poorly known, even in the 
relatively accessible southeastern part of the biome. 


INTRODUCTION 


Species of Neacomys Thomas, 1900, are small (<42 g) Neotropical rodents belonging to the 
tribe Oryzomyini in the cricetid subfamily Sigmodontinae. Commonly known as spiny mice or 
bristly mice, they occur in forested habitats of Central and South America, from eastern Panama 
to central Bolivia and central/western Brazil (Weksler and Bonyicino, 2015), Fourteen nominal 
taxa are referred to the genus, of which 12 are currently recognized as valid: N, amoenus Thomas, 
1903; N, dubosti Voss et al., 2001; N, guianae Thomas, 1905; N. macedoruizi Sanchez-Vendizu et 
al., 2018; N. minutus Patton et al., 2000; N. musseri Patton et al., 2000; N. paracou Voss et al., 2001; 
N. pictus Goldman, 1912; N. rosalindae Sanchez-Vendizt et al., 2018; N. spinosus (Thomas, 1882); 
N. tenuipes Thomas, 1900; and N. vargasllosai Hurtado and Pacheco, 2017. The nominal taxon 
pusillus Allen, 1912, has been considered a junior synonym of N. tenuipes (see Weksler and Bon- 
vicino, 2015), and the nominal taxon carceleni Hershkovitz, 1940, is currently considered to be a 
valid subspecies of N. amoenus (see Hurtado and Pacheco, 2017). 

In their monographic study of the small mammals from the Jurua River, Brazil, a right- 
bank tributary of the Amazon River, Patton et al, (2000) estimated the phylogenetic affinities 
of local populations of Neacomys and assessed their taxonomic status. Based on sequence data 
from one mitochondrial gene (cytochrome b), linear body measurements, and discrete mor- 
phological data, these authors revealed the existence of several previously unrecognized clades 
within the genus. In addition to providing the first gene tree for the group, these authors docu- 
mented high levels of genetic differentiation among several morphologically diagnosable lin- 
eages. Based on these results, Patton et al. (2000) described two new species from western 
Amazonia and suggested the existence of an undescribed taxon from eastern Amazonia, which 
they referred to as “Neacomys sp. clade 7." 

In a monograph on the mammals of Paracou, French Guiana, Voss et al. (2001) assessed 
the taxonomic status of local samples of Neacomys, provided an emended morphological diag- 
nosis for the genus, and described two new sympatric species for the Guiana Region. More 
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recently, three new species have been described from western Amazonia (eastern Peru) based 
on molecular, morphological, and karyotypic data (Hurtado and Pacheco, 2017; Sanchez- 
Vendizu et al., 2018). 

Cytogenetic studies have also indicated the existence of undescribed species of Neacomys 
in the eastern Amazon of Brazil. Silva et al. (2015) reported new karyotypes for specimens from 
Marajo Island and from Maraba, both in eastern Para state, and Oliveira da Silva et al. (2017) 
described two additional new karyotypes associated with specimens from the right and left 
banks of the Tapajés River. By using comparative chromosome painting analyses, Oliveira da 
Silva et al. (2019) described yet another new karyotype for specimens from Santa Barbara (in 
eastern Para) and suggested that this and the other cytotypes previously reported by Silva et 
al. (2015) and Oliveira da Silva et al. (2017) could represent new species. 

In view of the seven new species described in the past 20 years (Patton et al., 2000; Voss 
et al., 2001; Hurtado and Pacheco, 2017; and Sanchez-Vendizu et al., 2018), and the possible 
existence of undescribed taxa as suggested by molecular and cytogenetic studies (e.g., Patton 
et al., 2000; Hurtado and Pacheco, 2017; Sanchez-Vendizu et al., 2018; Oliveira da Silva et al., 
2019), it is clear that Neacomys is much more diverse than had been assumed by mid-20th 
century authors (e.g., Cabrera, 1961). However, despite the geographically restricted studies 
mentioned earlier, no comprehensive taxonomic revision has ever been undertaken for the 
genus. As a result, the systematics of Neacomys remains a work in progress, and little is yet 
known about the geographic limits of the species currently recognized as valid; additionally, 
phylogenetic relationships, especially of the unnamed forms in the eastern Amazon, remain 
obscure. In fact, the need for a critical taxonomic reevaluation of the genus is widely recognized 
(Woodman et al., 1991; Voss et al., 2001; Musser and Carleton, 2005; Catzeflis and Tilak, 2009; 
Weksler and Bonyicino, 2015; Hurtado and Pacheco, 2017; Sanchez-Vendizu et al., 2018). 

In this report, we review the alpha taxonomy of Neacomys from southeastern Amazonia 
and provide preliminary hypotheses about their phylogenetic relationships. Based on phyloge- 
netic analyses of the mitochondrial cytochrome b gene, together with external and craniodental 
characters, we define the species limits of the populations occurring in the southeastern parts 
of Brazilian Amazonia, In total, we recognize four species of Neacomys in southeastern Ama- 
zonia, three of them herein described as new. Our work represents a first step toward a thor- 
ough systematic review, offering novel data to enlighten future biogeographic, taxonomic, and 
evolutionary studies of this neglected group of rodents. 


Materials and Methods 


SOURCE OF MATERIAL: We analyzed a total of 174 specimens (skulls, skins, and entire 
fluid-preserved specimens), including name-bearing types and other important material from 
most of the distributional range of Neacomys (table 1; appendix 1). This examined material, as 
well as unexamined vouchers for specimens with analyzed molecular sequences, is deposited 
in the following collections (abbreviations in parentheses); American Museum of Natural His- 
tory, New York (AMNH); Natural History Museum, London (BMNH); Carnegie Museum, 
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TABLE 1. List of sequenced specimens of Neacomys included in phylogenetic analyses of cytochrome b. 
(Voucher material marked with asterisks (*) was examined by the authors.) 


Taxon 

amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus amoenus 
amoenus carceleni 
amoenus carceleni 
amoenus carceleni 
amoenus carceleni 
amoenus carceleni 
amoenus carceleni 
dubosti 

dubosti 

dubosti 


guianae 


Voucher 
UFMT 1669* 
UFMT 1659* 
UFMT 1367* 
UFMT 1374* 
UFMT 1373* 
UFMT 1370* 
MAS NB 46* 
MSA FM 02* 
INPA 3059* 
INPA 3064* 
INPA 3063* 
UFMT 3382* 
UFMT 3380* 
MVZ 155015* 
MVZ 190634* 
MVZ 190635* 
INPA 3057* 
MVZ 190372* 
USNM 588096 
USNM 588051 
USNM 584543 
USNM 584544 
INPA 3058* 
INPA 3060* 
INPA 3062* 
UFMT 1763* 
UEFMT 1757* 
ROM 104474 
ROM 105278 
ROM 105282 
ROM 105290 
ROM 105264 
USNM 574567 
CN 240* 

CM 76846 # 
V-1134 

CM 76847* 


GenBank 
MT462078 
MT462079 
MT462080 
MT462081 
MT462082 
MT462083 
MT462084 
MT462085 
MT462015 
MT462016 
MT462017 
MT462086 
MT462087 
KX792049 
MT462019 
MT462020 
MT462021 
KX792040 
KX792042 
KX792041 
KX792022 
MT462022 
KY859734 
KX792032 
KX792033 
MT462076 
MT462077 
KX792045 
KX792046 
KY859738 
KX792047 
MT462018 
KX792048 
MT462036 
FM210781 
FM210773 
FM210778 


Locality 

Brazil: Mato Grosso 
Brazil: Mato Grosso 
Brazil: Mato Grosso 
Brazil: Mato Grosso 
Brazil: Mato Grosso 
Brazil: Mato Grosso 
Brazil: Rondonia 
Brazil: Rondonia 
Brazil: Acre 

Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Mato Grosso 
Brazil: Mato Grosso 
Peru: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Peru: Cusco 

Peru: Cusco 
Bolivia: Santa Cruz 
Bolivia: Santa Cruz 
Brazil: Acre 

Brazil: Acre 

Brazil: Acre 

Brazil: Mato Grosso 
Brazil: Mato Grosso 
Ecuador: Napo 
Ecuador: Napo 
Ecuador: Napo 
Ecuador: Napo 
Ecuador: Napo 
Ecuador: Pastaza 
Brazil: Para 
Suriname: Nickerie 
France: French Guiana 


Suriname: Nickeri 


Source 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

Hurtado and Pacheco, 2017 
Patton et al., 2000 * 

Patton et al., 2000 * 

This study 

Patton et al., 2000 

Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Patton et al., 2000 

Patton et al., 2000 * 
Sanchez-Vendizt et al., 2018 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
This study 

This study 

Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Sanchez-Vendizt et al., 2018 
Hurtado and Pacheco, 2017 
Patton et al., 2000 * 
Hurtado and Pacheco, 2017 
This study 

Catzeflis et al., 2009 
Catzeflis et al., 2009 
Catzeflis et al., 2009 


2020 


Taxon 
guianae 
guianae 
macedoruizi 
macedoruizi 
marajoara 
marajoara 
marajoara 
marajoara 
marajoara 
marajoara 
marajoara 
marajoara 


marajoara 


minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 
minutus “downriver” 


minutus “downriver” 


minutus “upriver” 
minutus “upriver” 
minutus “upriver” 
minutus “upriver” 
minutus “upriver” 
minutus “upriver” 
minutus “upriver” 
musseri 

mussert 


musseri 


Voucher 

INPA 7102* 
INPA 7100* 
MUSM 45054? 


MUSM 45053 * 


MPEG40440 
MPEG 40443 
MPEG 40441 
MPEG 40446 
MPEG 40435 
MPEG 40434 
MPEG 404324 
MPEG 40431 
MPEG 40429 
SISPUR 159* 
SISPUR 151* 
INPA 5392* 
INPA 5387* 
INPA 3048* 
INPA 3049* 
EE 110* 
MVZ 190360 
MVZ 190361 
MVZ 191209 
MVZ 190359 
MVZ 190363 
INPA 3051* 
INPA 3047* 
INPA 2689* © 
INPA 3055* 
INPA 3053* 
MVZ 190358 
INPA 3050* 
MVZ 190362 
INPA 3891* 
INPA 3056* 
MVZ 171487 
MVZ 171488 
INPA 3046* ! 


GenBank 
MT462037 
MT462038 
KY859731 
KY859732 
KX752075 
KX752072 
MT462067 
KX752080 
MT462068 
MT462069 
MT462070 
MT462071 
MT462072 
MT462039 
MT462040 
MT462023 
MT462007 
MT462008 
MT462009 
MT462041 
KY859739 
KX792069 
KX792072 
KX792068 
KX792071 
KX792065 
KX792063 
MT462010 
MT462011 
MT462012 
U58392 
KX792064 
KX792070 
KX792066 
KY886326 
KX792074 
KY859742 
U58392 


Locality 

Brazil: Amazonas 
Brazil: Amazonas 
Peru: Huanuco 
Peru: Huanuco 
Brazil: Para 
Brazil: Para 
Brazil; Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil; Para 
Brazil: Para 
Brazil: Para 
Brazil: Amazonas 
Brazil; Amazonas 
Brazil: Amazonas 
Brazil; Amazonas 
Brazil: Amazonas 
Brazil; Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Brazil: Amazonas 
Peru: Cusco 

Peru: Cusco 


Brazil: Acre 
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Source 
This study 
This study 


Sanchez-Vendizt et al., 2018 
Sanchez-Vendizu et al., 2018 
Oliveira da Silva et al., 2017 
Oliveira da Silva et al., 2017 


This study 


Oliveira da Silva et al., 2017 


This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
Patton et al., 2000 


Patton et al., 2000 
Patton et al., 2000 
Patton et al., 2000 
Patton et al., 2000 


Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 


This study 
This study 
This study 


Sanchez-Vendizu et al., 2018 
Sanchez-Vendizt et al., 2018 
Sanchez-Vendizu et al., 2018 
Sanchez-Vendizu et al., 2018 
Sanchez-Vendizu et al., 2018 
Hurtado and Pacheco, 2017 


Patton et al., 2000 


Hurtado and Pacheco, 2017 


Taxon 
paracou 
paracou 
paracou 
paracou 
paracou 
paracou 
paracou 
paracou 
paracou 
paracou 
paracou 
rosalindae 
rosalindae 
rosalindae 
rosalindae 
rosalindae 
rosalindae 
rosalindae 
species 1 
species | 
species 1 
species | 
species 1 
species 2 
species 2 
species 2 
species 2 
species 2 
species 2 
species 2 
species 2 
spinosus 


tenuipes 


vargasllosai 
vargasllosai 
vargasllosai 
vossi 


vosst 
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Voucher 

CN 279* 

CN 263* 
MPEG 40080* 
MPEG 40078* 
CN 129* 
MPEG 40000* 
MPEG 39998* 
INPA 7089 * 
INPA 7138 * 
ROM 114325* 
ROM 114315* 
ROM 104560 
ROM 105265 
ROM 105314 
ROM 105315 
MVZ 153530 
MVZ 155299 
KU 158172 
INPA 4190 
INPA 4191 
INPA 4189 
INPA 4192 
MPEG 45483 
MPEG 42838 
JUR 82 

JUR 19 

UFPA 1227 
UFPA 1530 
UFPA 1413 
JUR 042 
MPEG 42901 
MUSM 36924 


BMNH 
1899.10.3,34*8 


MVZ 172650 
MVZ 172654 
MVZ 172655 
UFPA 1277 
UFPA 1284 


GenBank 
MT462042 
MT462043 
MT462044 
MT462045 
MT462046 
MT462047 
MT462048 
MT462049 
MT462050 
FM210782 
FM210767 
KX792050 
KX792051 
KX792052 
KX792053 
KX792054 
KY859730 
KX792055 
KX792061 
KX792062 
KX792059 
KX792060 
MT462051 
MT462052 
MT462053 
MT462054 
MT462034 
MT462035 
MT462055 
MT462056 
MT462057 
KX258228 
KX792081 


KX792082 
MT462013 
MT462014 
MT462024 
MT462025 


Locality 


Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 


Brazil: 


Para 
Para 
Para 
Para 
Para 
Para 
Para 
Amazonas 


Amazonas 


Surinam: Brokopondo 


Surinam; Brokopondo 


Ecuador: Napo 


Ecuador; Napo 


Ecuador; Napo 


Ecuador: Napo 


Peru: Amazonas 


Peru: Amazonas 


Peru; Loreto 


Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 
Brazil: 


Amazonas 
Amazonas 
Amazonas 
Amazonas 
Amazonas 
Para 
Para 
Para 
Para 
Para 
Para 
Para 


Para 


Peru: Amazonas 


Colombia: Cundina- 


Marca 


Peru: Puno 
Peru; Puno 
Peru: Puno 
Brazil; Para 


Brazil: Para 
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Source 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

Catzeflis et al., 2009 
Catzeflis et al., 2009 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Sanchez-Vendizt et al., 2018 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
Hurtado and Pacheco, 2017 
This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

Hurtado and Pacheco, 2017 
Patton et al., 2000 


Hurtado and Pacheco, 2017 
Patton et al., 2000 * 

Patton et al., 2000 * 

This study 

This study 


2020 SEMEDO ET AL.: SYSTEMATICS OF NEACOMYS THOMAS, 1900 


Taxon Voucher 
vossi UFPA 1647 
vossi UFPA 1467 
vossi UFPA 1583 
vossi UFPA 1736 
vossi UFPA 1577 
vossi UFPA 1691 
vossi UFPA 1520 
vossi UFPA 1444 
vossi UFPA 1391 
vossi UFPA 1654 
vossi UEFPA 1487 
xingu UFMT 1275 
xingu UFMT 1273 
xingu UFMT 1268! 
xingu MPEG 41805 
xingu MPEG 42019 
xingu MPEG 41804 
xingu MPEG 41996 
xingu PSA 69 
xingu MPEG 41991 
xingu PSA 46 

xingu USNM 549553 } 


* Paratype of Neacomys dubosti. 

> Paratype of Neacomys macedoruizi. 
“Holotype of Neacomys macedoruizi. 

4 Holotype of Neacomys marajoara. 

© Holotype of Neacomys minutus. 

* Paratype of Neacomys musseri. 

8 Holotype of Neacomys tenuipes. 

h Holotype of Neacomys vossi. 

' Holotype of Neacomys xingu. 

) Neacomys “clade 7” of Patton et al., 2000. 


« Unpublished data from Patton et al., 2000. 


GenBank 

MT462026 
MT462027 
MT462028 
MT462073 
MT462029 
MT462030 
MT462031 
MT462074 
MT462075 
MT462032 
MT462033 
MT462058 
MT462059 
MT462060 
MT462061 
MT462062 
KX752073 
MT462063 
MT462064 
MT462065 
MT462066 
KX792080 


Locality 

Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 
Brazil: Para 


Brazil: Para 


Source 

This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 


Oliveira da Silva et al., 2017 


This study 
This study 
This study 
This study 


Hurtado and Pacheco, 2017 


oa 
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Pittsburgh (CM); Instituto Nacional de Pesquisas da Amaz6nia, Manaus (INPA); Pontificia 
Universidade Catdélica de Minas Gerais, Belo Horizonte (MCN-M); University of Kansas Bio- 
diversity Research Center, Lawrence (KU); Muséum d Histoire Naturelle de la Ville de Geneve, 
Geneva (MHNG), Muséum National d'Histoire Naturelle, Paris (MNHN); Museu Paraense 
Emilio Goeldi, Belém (MPEG); Museo de Historia Natural de la Universidad Nacional Mayor 
de San Marcos, Lima (MUSM); Museum of Vertebrate Zoology, University of California, 
Berkeley (MVZ); Royal Ontario Museum, Toronto (ROM); Colec4éo Zooldgica da Universidade 
Federal de Mato Grosso, Cuiaba (UFMT); Universidade Federal do Para, Belém (UFPA); and 
National Museum of Natural History, Smithsonian Institution, Washington (USNM). Speci- 
mens with field numbers prefixed by CN, JUR, and PSA are currently held at the MPEG and 
by MAS NB, MSA FM, EE, and SISPUR are currently held at INPA. Those specimens will be 
deposited in the aforementioned collections. 

CoLLecTING Locatitizs: Geographic data for collection localities were obtained from 
specimen tags, For specimens unaccompanied by geographic coordinates, we referred to Payn- 
ter and Traylor (1991), Paynter (1993), Prado and Percequillo (2013), Gardner (2008), and 
Patton et al. (2015) to obtain this information. 


Molecular Analyses 


TAXON SAMPLING: We obtained partial sequences (801 bp) of the mitochondrial cytochrome b 
gene (Cytb) trom 75 Amazonian specimens of Neacomys. These new sequences are deposited in 
GenBank with accession numbers MT462007-MT462087, In addition, we downloaded 54 ingroup 
(Neacomys) sequences from GenBank and another six unpublished sequences from Patton et alls 
(2000) dataset (table 1). In total, our ingroup dataset includes all currently recognized species in the 
genus, with the sole exception of N. pictus (for which no sequence data are available), Among other 
nomenclaturally important sequences, we analyzed Cytb data from the holotypes of N. minutus 
(INPA 2689) and N. tenuipes (BMNH 1899,10,3.34), a paratype of N. musseri (INPA 3046), and the 
holotype (MUSM 45053) and a paratype (MUSM 45054) of N. macedoruizi. 

We obtained sequences of the following outgroup taxa from Genbank (accession numbers 
in parentheses): Thomasomys baeops (DQ914654), Hylaeamys megacephalus (KF815441), Oli- 
goryzomys microtis (OMU58381), Scolomys ucayalensis (EU579518), Oreoryzomys balneator 
(EU579510), and Microryzomys minutus (AF108698). The T. baeops sequence was used to root 
the trees in our phylogenetic analyses. Outgroup taxa were chosen based on the phylogenies 
of Weksler (2003), Weksler (2006), and Percequillo et al. (2011). 

LABORATORY PROCEDURES: We extracted genomic DNA from ethanol-preserved tissues 
using the salt-extraction method (Aljanabi and Martinez, 1997). The Cytb gene was PCR-ampli- 
fied using primers developed by Smith and Patton (1993): MVZ05 (5'-CGAAGCTTGATAT- 
GAAAAACCATCGTTG-3’) and MVZ16 (5’AAATAGGAARTATCAYTCTGGTTTRAT-3’), 
Each PCR reaction contained 2.5 LL of 10x buffer, 1.0 pL of MgCl, (50 mM), 1.25 pL of dNTP 
mix (1.25 mM for each nucleotide), 1 uL of each primer (10 mM), 0.2 L of Ludwig Biotech- 
nolgia Taq DNA Polymerase (5 p/uL), 1 uL of DNA template, and water to a final volume of 


2020 SEMEDO ET AL.; SYSTEMATICS OF NEACOMYS THOMAS, 1900 9 


25 tL. For some difficult templates, concentrations of reagents were adjusted to increase PCR 
success. PCR conditions used a preheating step of initial denaturation at 94° C for 3-5 min, 
followed by 35 cycles at 94° C for 30 sec, 49° C for 45 sec, and 72° C for 1 min and 45 sec, 
and a final extension at 72° C for 5-7 min. PCR products were cleaned using Exonuclease | 
and shrimp alkaline phosphatase (Fermentas), and sequenced using an ABI 3130xl Genetic 
Analyzer (Life Technologies), 

MOLECULAR Data ANALYSES: We aligned sequences with Clustal W (Thompson et al., 
1994) using default parameters. Phylogenetic analyses were performed using maximum likeli- 
hood (ML) and Bayesian inference (BI). We coded all missing bases as unknown (“?”) prior to 
phylogenetic analyses. PartitionFinder ver. 2.1.1 (Lanfear et al., 2017)—which heavily relies on 
PhyML (Guindon et al., 2010)—was employed to select both the most suitable partition scheme 
and the best-fit model of nucleotide substitution for each data subset, The Akaike Information 
Criterion (AIC) was used to select the best-fit models; only models that could be applied in 
MrBayes were considered, and the “greedy” algorithm was implemented. 

In order to find the best topology according to the ML optimality criterion, we conducted 
an analysis with 50 independent searches in GARLI 2.01 (Zwickl, 2006) using default settings. 
Maximum-likelihood bootstrap analysis was also carried out in GARLI 2.01 and implemented 
in the CIPRES Science Gateway (Miller et al., 2010) using 100 pseudoreplicated data matrices 
with 10 searches performed on each. The Bayesian analysis was conducted using the Markov 
Chain Monte Carlo (MCMC) sampling approach in MrBayes v. 3.2.3 (Huelsenbeck and Ron- 
quist, 2001; Ronquist and Huelsenbeck, 2003; Altekar et al., 2004; Ronquist et al., 2012) and 
was also implemented through the CIPRES Science Gateway. The search started with a random 
tree and consisted of one cold chain and three heated chains. The Markov chains were run for 
100,000,000 generations, and trees were sampled every 1000 generations. Default values were 
kept for the “relburnin” and “burninfrac” options in MrBayes; hence, the first 25,000,000 gen- 
erations (i.e., 250,000 trees or 25%) were discarded as burn-in, and posterior probability esti- 
mates of all model parameters were based on the remaining (750,000) trees. 

We calculated average pairwise genetic distances within and between groups using the 
uncorrected (p-distance) option in MEGA 6 (Tamura et al., 2011). Uncorrected distances were 
chosen to allow comparison with divergence estimates reported in previous studies of Neaco- 
mys (e.g., by Hurtado and Pacheco, 2017; Sanchez-Vendizu et al., 2018), 


Morphological and Morphometric Analyses 


MorpxHo.ocy: We examined the following external and craniodental characters for taxo- 
nomic variation: dorsal and ventral pelage color, tail color, morphology of the palate, carotid 
circulation type, features of the auditory bullae, and dental morphology. We made comparisons 
predominantly among specimens of the same age class to avoid conflating ontogenetic variation 
with taxonomic differences. Nomenclature for external and craniodental morphology follows 
Voss (1988), Patton et al. (2000), Voss et al. (2001), Weksler (2006), Hurtado and Pacheco 
(2017), and Sanchez-Vendizu et al. (2018). 
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AGE Criteria: In order to distinguish ontogenetic from taxonomic variation, we grouped 
specimens into age classes based on patterns of molar eruption and degree of tooth wear (fig. 
1), using the criteria defined by Voss (1991): 


Age class 1 M3 incompletely erupted or unworn; 

Age class 2 M3 fully erupted and exhibiting slight to moderate wear (some dentine exposed), 
but the occlusal surface still tubercular (the paracone raised and prominent), not flat; 
Age class 3. M3 well worn, the occlusal surface flat or concave; M1-2 tubercular (the major 

cusps all separate and prominent); anteroloph of M2 distinct, not fused with paracone; 
Age class 4 M3 flat or concave; M1-2 with cusps worn almost or quite flat but not below 
widest part of crown; anteroloph of M2 obliterated, fused with paracone; 
Age class5 M1-3 all worn flat or concave, below widest part of crowns; most details of occlu- 
sal topography obliterated. 


Hereafter, we refer to specimens in age class 1 as “juveniles,” specimens in age classes 2-4 
as “adults,” and specimens in age class 5 as “old adults.” 

MEASUREMENTS: We obtained 21 craniodental measurements (in millimeters, mm) that 
have been previously used in taxonomic studies of Neotropical cricetids (Voss, 1988; Patton et 
al., 2000; Voss et al., 2001; Flores et al., 2010). The following dimensions were measured with 
digital calipers to the nearest 0.01 mm while skulls were examined at low magnification under 
a stereomicroscope: 


BH _ Braincase height, measured in the midline from the basisphenoid-basioccipital suture 
on the ventral surface of the skull to the frontal-parietal suture on the dorsal surface; 

BI Breadth of incisive foramen, the greatest transverse dimension across both incisive 
foramina; 

BM1_ Breadth of M1, greatest crown breadth of the first upper molar (M1); 

BPR Breadth of palate at rostrum, measured between most posterior lower edges of the left 
and right infraorbital foramina on the ventral side of skull; 

BZP Breadth of zygomatic plate, least distance between the anterior and posterior edges of 
the zygomatic plate; 

BB Bullar breadth, distance between anterior opening of carotid foramen to ectotimpanic 
dorsal process; 

CIL Condylo-incisive length, from the greater curvature of one upper incisor to the articu- 
lar surface of the condyle on the same side; 

LCIB Least condyloid-incisor breadth, greatest distance from the lower incisor base to the 
posterior margin of the condyloid process on the same side of the mandible; 

LIB _ Least interorbital breadth, least distance across the frontal bones between the orbital 
fossae; 

LD Length of diastema, from the crown of M1 to the lesser curvature of the incisor on the 
same side; 

LIF Length of incisive foramen, greatest anterior-posterior dimension of one incisive 
foramen; 
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FIG. 1. SEM images of molar toothwear age classes of Neacomys amoenus used in this study (based on Voss, 
1991). From left to right: age class 1 (UFMT 4050), age class 2 (UFMT 1758), age class 3 (UFMT 1760), age 
class 4 (UFMT 1759), and age class 5 (UFMT 915). Arrows indicate the incompletely erupted M3 in age class 
1, the distinct M2 anteroloph (not fused with the paracone) in age class 3, and the obliterated M2 anteroloph 
(fused with the paracone) in age class 4. 


LLD Length of lower diastema, from the crown of m1 to the lesser curvature of the incisor 
on the same side of the mandible; 

LLM Length of lower molars, crown length from m1 to m3; 

LM Length of molars, crown length from M1 to M3; 

LN Length of nasals, greatest anterior-posterior dimension of one nasal bone; 

LPB Length of palatal bridge, from the posterior margin of the upper incisors to the ante- 
rior margin of the mesopterygoid fossa; 

MB __Mastoid breadth, distance across the cranium at the mastoid processes; 

OL —_ Orbital length, internal distance between the anterior and posterior margins of the right 
or the left orbit; 


RB Rostral breadth, distance across the outside margins of the left and right nasolacrimal 
capsule; 

RL Rostral length, diagonal measurement taken from the anterior margin of the orbit to 
the anterior margin of the nasal bone on the same side; 

ZB Zygomatic breadth, greatest transverse dimension across the squamosal zygomatic 
processes. 


In addition to these craniodental measurements, we transcribed five external measure- 
ments (in mm) taken in the field by collectors and recorded on specimen labels: total length 
(TL); length of head and body (HBL); length of tail (LT); length of hind foot (HF); and length 
of ear (Ear). We also transcribed weight (in grams, g) as recorded by collectors. 

STATISTICAL ANALYSES: We calculated standard descriptive univariate statistics for all exter- 
nal and craniodental measurements. In order to minimize potentially confounding ontogenetic 
variation, we included only adult and old adult specimens (age classes 2-5) in these analyses. 
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Principal component analyses (PCA) were conducted to assess whether morphometric data 
are congruent with results from the molecular analyses, qualitative morphological compari- 
sons, and previously published karyotypic studies. For PCA, measurements were log-trans- 
formed and principal components were extracted from the variance-covariance matrix. 
Statistical analyses were performed using SPSS 22.0 for Windows. 


RESULTS 
Molecular Analyses 


The Cytb gene fragment we sequenced was 801 base pairs long and included 470 conserved 
sites, 331 variable sites, and 287 parsimony-informative sites. Monophyly of Neacomys was 
recovered with strong support from our Bayesian analysis (Bayesian posterior probability 
[BPP] = 0.97); whereas our ML analysis provided only marginal support (bootstrap support 
[BS] = 51%) for generic monophyly (fig. 2A-C). 

The gene trees obtained from our Bayesian and ML analyses both recovered 17 distinct 
lineages of Neacomys, 11 of which could be associated with taxa that are currently recognized 
as valid species: N. dubosti, N. guianae, N. minutus, N. musseri, N. paracou, N. spinosus, N. 
rosalindae, N. macedoruizi, N. amoenus, N. vargasllosai, and N. tenuipes. Among the remaining 
six lineages, three are formally described in this report (N. vossi, N. xingu, and N. marajoara) 
and three others are referred to by informal nomenclature: N. minutus “upriver” and N. minu- 
tus “downriver” (after Patton et al., 2000), Neacomys “sp. 1,’ and Neacomys “sp. 2.” All these 
lineages received BPP support values > 0.96, except for N. minutus “downriver” (BPP = 0.82). 
All undescribed lineages also had bootstrap values = 94%, except for N, xingu and N. marajoara 
(with bootstrap values of 58% and 62%, respectively; fig. 2B). Mean interspecific genetic dis- 
tances ranged from 3.9% (between N. marajoara and N. xingu) to 15.58% (between N. paracou 
and N. musseri; table 2), whereas intragroup divergences ranged from 0% in N. macedoruizi 
(two sequences from one locality) and N. vargasllosai (three sequences from one locality) to 
1.88% in N. amoenus (34 sequences from 16 localities). 

The only relationship that differs between our BI and ML topologies is the position of N. 
tenuipes, which is sister to a negligibly supported clade containing N. minutus, N. macedoruizi, 
N. rosalindae, N. musseri, Neacomys sp. 1, and N. guianae in the BI tree. By contrast, N, tenuipes 
is sister to a negligibly supported clade containing N. dubosti, N. marajoara, Neacomys sp. 2, 
N. vossi, and N. xingu in the ML tree (appendix 2). 

Based on the BI topology, and for the purposes of this report, we recognize four groups of 
species within Neacomys. Two of these groups correspond to the Paracou and Spinosus groups 
(fig. 2C) of Hurtado and Pacheco (2017), whereas the Tenuipes Group (fig. 2A) of Hurtado and 
Pacheco (2017) was recovered as two monophyletic groups in our analysis. Herein we restrict the 
Tenuipes Group to N, tenuipes, N, guianae, Neacomys sp. 1, N, musseri, N, rosalindae, N, mace- 
doruizi, and N, minutus; additionally, we recognize a Dubosti Group for N. dubosti, Neacomys sp. 
2, and the three new species described in this report (N. marajoara, N. xingu, and N. vossi). 
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Neacomys paracou (comprising the monotypic Paracou Group) was recovered as sister to the 
remaining species of Neacomys with strong Bayesian support, whereas the Spinosus Group was 
recovered as sister to the Dubosti and Tenuipes groups with strong Bayesian support. The sister- 
group relationship between the Dubosti and Tenuipes groups also received strong Bayesian support. 
Within the Spinosus Group, N. vargasllosai was recovered as sister to N. spinosus + N, amoenus with 
consistently strong support. Within the Dubosti Group, the sister-group relationship of N. dubosti 
to other group members in the nested clades (N. marajoara (N. xingu (Neacomys sp. 2 + N, vossi))) 
is also strongly supported. Lastly, with negligible support values, the Tenuipes Group includes N. 
tenuipes as sister to ((N. guianae + Neacomys sp. 1) ((N. musseri + N. rosalindae) ((N. macedoruizi 
+ N. minutus “upriver”) N. minutus “downriver’))). In general, Cytb provided satisfactory resolution 
of supraspecific phylogenetic relationships. However, among the new species described herein, the 
relationships of N, xingu and N. vossi had negligible support values. 


Morphometric Analyses 


Descriptive statistics for adult and old adult specimens of Neacomys dubosti and three new 
species from southeastern Amazonia are shown in table 3. In general, our results document 
univariate and multivariate morphometric overlap among all species in the Dubosti Group. In 
all pairwise multivariate analyses herein provided, the first principal components comprise 
only positive coefficient values, indicating that they are predominantly related to variation in 
skull size. Comparing N. dubosti with N. marajoara, the plot of first (PC1) and second principal 
components (PC2) shows that specimens of N, dubosti group at the higher extremity of the size 
axis (fig. 3A, appendix 3), The PC1 corresponds to 46.98% of the total variation and is strongly 
related to CIL and ZB while the PC2 corresponds to 12.62% of the total variation and is 
strongly related to LIF and LM. Similarly, N. dubosti tended to occupy the higher extremity of 
the size axis when compared to N. vossi (fig. 3B, appendix 3), In this analysis, PC1 corresponds 
to 47,22 % of the total variation and is strongly related to CIL and ZB, while the PC2 corre- 
sponds to 11.31% of the total variation and is strongly related to LIF and LCIB. By contrast, N. 
dubosti largely overlaps N. xingu in the morphospace (fig. 3C, appendix 3). In this analysis, the 
PC1 explains 35.58% of the total variation and is strongly related to CIL, LD and ZB, while the 
PC2 explains 13.33% and is strongly related to BH, LN, and LCIB. 

When N. marajoara was compared with N. vossi (fig. 3D, appendix 3), the PC1 explained 28.97 
% of the total variation, while the PC2 explained 15.48%. The specimens largely overlapped on both 
axes, and the measurements that most contributed for the variation were CIL, LIF, LN, RL, and ZB. 
By comparing N. marajoara and N. xingu (fig. 3E, appendix 3), the PC1 explained 41.20 % of the 
total variance, while the PC2 explained 12.35%. The latter species appeared at the higher extremity 
of the size axis, and the measurements with highest loadings were BZP, CIL, LN, RL, and ZB. 

In the PCA including N, vossi and N. xingu (fig. 3F, appendix 3), the PC1 explained 41,13 
% of the total variance, while the PC2 explained 10.94%, The latter species appeared at the 
higher extremity of the size axis, and the variables that better contributed for variation were 
CIL, LM, LN, RL, and ZB. 
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FIG. 2 A-C (above and following two pages). Bayesian phylogeny of Neacomys based on cytochrome b sequence 
data. Numbers represent Bayesian posterior probabilities (above) and ML bootstrap support (below) at nodes 
recovered by both analyses (bootstrap support for nodes only recovered by our ML analysis are shown in 
appendix 2). Only ML bootstrap support >50% is shown. Branch tips are labeled with geographic identifiers 
and voucher codes listed in table 1. Specimens marked with asterisks (*) indicate holotypes. 
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TABLE 3. Descriptive statistics for external and craniodental dimensions (in mm) and mass (in g) of species in the 


Dubosti Group of Neacomys. 


one 


HBL 


LT 


HEF 


Ear 


Mass 


BH 


BIF 


BM1 


BPR 


BZP 


BB 


CIL 


LCIB 


LIB 


LD 


LIF 


LLD 


LLM 


LM 


LN 


dubosti 

157 +7 
148-167 (9) 
79+4 
74-88 (9) 
77+A4 
71-83 (9) 
21+1 
20-23 (9) 
eo 
12-15 (9) 
15+3 
12-20 (9) 
714+ 0.27 
6.61-7.48 (12) 


1.51+0.1 
1,39-1.68 (12) 


0.90 + 0.04 
0.82-0.98 (12) 


4.02 + 0.18 
3.73-4,32 (12) 


2.06 + 0.19 
1.86-2.58 (12) 


3.66 + 0.18 
3,38-3.93 (12) 


18.57 + 0.59 
17.38-19.43 (12) 


8.13 + 0.97 
5,33-8.96 (12) 


4.65 + 0.12 
4.47-4.87 (12) 


5.36 + 0.25 
4,94-5.69 (12) 


2.98 + 0.26 
2.64-3.47 (12) 


3.0 + 0.30 
2.24-3.3 (12) 


2.97 + 0.13 
2.56-3.06 (12) 


2.75 + 0.07 
2.63-2.84 (12) 


7.81 + 0.35 
7.23-8.51 (12) 


marajoara 
144 + 10 
128-153 (9) 


72+4 
62-77 (9) 
72 +8 
55-79 (9) 
17 +2 
13-19 (9) 
12 +1 
11-14 (9) 
1Z*.2 
9-16.5 (9) 
6.7 + 0.32 
6.33-7.08 (7) 


1.40 + 0.08 
1.30-1.55 (7) 


0.85 + 0.03 
0.81-0.88 (7) 


3.76 + 0.27 
3.37-4.06 (7) 


1,84 + 0.11 
1.67-1.97 (7) 


3.46 + 0.22 
3.17-3.8 (7) 
17.51 + 0.36 
16.92-18.06 (7) 


7.71 40.3 
7.13-7.97 (7) 


4,47 + 0.2 
4,23-4.74 (7) 


5.01 + 0.20 
4.69-5.28 (7) 


2.99 + 0.10 
2.83-3.16 (7) 


2.89 + 0.23 
2.47-3.14 (7) 


2.77 + 0.04 
2.71-2.82 (7) 


2.53 + 0.06 
2.42-2.61 (7) 


7.45 + 0.45 
6.62-7.90 (7) 


xingu 
146 +5 
138-156 (7) 


72 
68-75 (7) 
7445 
68-85 (7) 

20 + 3 

14-26 (7) 
2+1 
11,5-14.5 (6) 


1443 
11.5-18.5 (7) 


6.77 + 0.2 
6.4—7 (13) 
1.41 + 0.11 
1.21-1.58 (13) 


0.87 + 0.04 
0.80-1.01 (13) 


4.09 + 0.20 
3.70-4.47 (13) 


1,99 + 0.23 
1.73-2.62 (13) 


3.63 + 0.17 
3,39-3.92 (13) 


18.06 + 0.58 
17.44-19.28 (13) 


8.17 + 0.53 
7.34-9.36 (13) 


4.56 + 0.26 
4.11-4.93 (13) 


5.06 + 0.15 
4.74-5.37 (13) 


2.98 + 0.23 
2.48-3.26 (13) 


2.87 + 0.18 
2.43-3.08 (13) 


2.92 + 0.08 
2.8-3.13 (13) 


2.75 + 0.09 
2.60-2.90 (13) 


7.89 + 0.46 
7.02-8.90 (13) 


Vosst 
140 + 13 
110-156 (14) 


68 + 10 
65-77 (14) 
72 +6 
60-84 (14) 
16 +3 
10-20 (14) 
12+3 
8-19 (14) 
12+2 
9-15.5 (14) 
6.71 + 0.38 
6.02-7.24 (8) 


1.44 + 0.09 
1.31-1.56 (8) 


0.85 + 0.08 
0.75-0.97 (8) 


3.92 + 0.27 
3,50-4.34 (8) 


1,88 + 0.07 
1.80-2.01 (8) 


3.51 + 0.11 
3.43-3.75 (8) 


17.37 + 0.48 
16.75-18.3 (8) 


7.93 + 0.27 
7.58-8.34 (8) 


4.46 + 0.18 
4.25-4.77 (8) 


4.92 + 0.14 
4,80-5.20 (8) 


2.91 + 0.30 
2.53-3.57 (8) 


3.34 + 0.99 
2.64-4.96 (8) 


2.77 + 0.08 
2.64-2.87 (8) 


2.58 + 0.08 
2.45-2.65 (8) 


7.24 + 0.39 
6.56-7.84 (8) 
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dubosti marajoara xingu vossi 
LPB 3.95 + 0,27 3.65 + 0.18 3.87 + 0.24 3.66 + 0,18 
3.49-4.45 (12) 3.31-3.85 (7) 3.40-4.25 (13) 3.46-4.06 (8) 
MB 8.86 + 0.23 8.27 + 0.22 8.51 + 0.22 8.34 + 0.12 
8.53-9.27 (12) 8.02-8.64 (7) 8.27-8.92 (13) 8.20-8.50 (8) 
OL 7.19 + 0.33 6.94 + 0.16 7.04 + 0.28 6.7 + 0.33 
6.45-7.76 (12) 6.74-7.12 (7) 6.51-7,44 (13) 6.16-7.28 (8) 
RB 3.77 + 0.24 3.80 + 0.20 3.89 + 0.30 3.77 + 0.16 
3.41-4.40 (12) 3.43-4,06 (7) 3.3-4.43 (13) 3.48-3.95 (8) 
RL 6.58 + 0.26 6.36 + 0.32 6.58 + 0.33 6.06 + 0.48 
6.12-6.98 (12) 5.84-6.68 (7) 6.09-7.20 (13) 5.43-6.82 (8) 
ZB 11.15 + 0.43 10.51 + 0.20 10.91 + 0.41 10.52 + 0.30 
10.46-11.78 (12) 10.21-10.76 (7) 9.93-11.51 (13) 10.13-11.05 (8) 
TAXONOMIC ACCOUNTS 


Our molecular analyses revealed the existence of 17 lineages of Neacomys with high levels 
of genetic divergence; of these, four occur in southeastern Amazonia (south of the Amazon 
River and east of the Rio Madeira). Although not strongly morphometrically differentiated, 
these previously unnamed southeastern lineages can be recognized as species based on our 
molecular analyses (which provide evidence of mtDNA divergence), by the karyotypic traits 
previously reported in the literature (which suggest they are reproductively isolated), and by 
qualitative morphological characters (which provide evidence of nuclear-gene divergence). 
Three of these species are described below, whereas the fourth (Neacomys “sp. 2”) will be 
described in another report with different authorship. Because all the new species belong to 
the Dubosti Group, our comparisons are restricted to members of that clade. Relevant sum- 
maries of morphometric variation and diagnostic traits are provided in tables 3 and 4, 
respectively. 


Neacomys marajoara, new species 


Marajoara Spiny Mouse 
Figures 4, 7 


Ho.ortype: The holotype (MPEG 40432) is an adult male (age class 3), collected on 19 
January 2009, by R.V. Rossi (field number MAJ 23; fig. 4) in a pitfall trap. The specimen consists 
of a stuffed skin (missing the tip of the tail), skull, and skeleton; a tissue sample of this speci- 
men is preserved in ethanol, and a partial cytochrome b sequence that we obtained from it has 
been deposited in Genbank with accession number MT462070. 

Type Locatity: Tauari Farm, municipality of Chaves, Marajo Island, state of Para, Brazil 
(0°39'S, 50°11'W, figs. 5, 6). 

Diacnosis: Neacomys marajoara is a small species (table 3) that differs from congeneric 
taxa by the following combination of craniodental traits: skull delicate; interorbital region rela- 
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FIG. 3. Scatter plots of PC1 (horizontal axis) and PC2 (vertical axis) for pairwise analyses of craniodental 
measurement data for species in the Dubosti Group of Neacomys. In each plot, species are represented by 
minimum convex polygons that enclose all points in the plane of PC1 and PC2. A, N. dubosti (dub) versus 
N. marajoara (mar); B, N. dubosti versus N. vossi (vos); C, N. dubosti versus N. xingu (xin); D, N. marajoara 
versus N. vossi; E, N. marajoara versus N. xingu; F, N. vossi versus N. xingu. 
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FIG. 4. The holotype of Neacomys marajoara, (MPEG 40432) in life. Photo by Rogério V. Rossi. 


tively broad; nasal bones straight anteriorly; supraorbital margins slightly convergent anteri- 
orly; subsquamosal fenestra large (almost half the size of the postglenoid foramen on each side 
of the skull); paraoccipital processes separated from the auditory bullae; sphenopalatine fora- 
men large; carotid circulation usually pattern 2 (sensu Voss, 1988); maxillary part of incisive 
septum (between the incisive foramina) narrow; M1 usually with flat and undivided antero- 
cone; and M1 anteroloph usually fused with anterolabial conule. 

MORPHOLOGICAL DESCRIPTION: Dorsal pelage dark brown finely sprinkled with orange 
(fig. 4); ventral pelage varying from pure white to yellowish white, separated from the dorsal 
pelage by a very thin orange lateral line. Superciliary, genal, and mystacial vibrissae blackish 
and long (extending behind ears when laid back alongside the head); submental vibrissae 
absent; interramal vibrissae short and white. Ears small and rounded; postauricular hairs gray 
based with orange tips, forming an orange tuft behind each pinna. Ungual tufts white, longer 
than claws in most specimens examined (except MPEG 40435 and MPEG 40446, in which 
ungual tufts are as long as the claws); fore- and hind feet covered dorsally with buffy-cream 
hairs; hind feet narrow and elongate with small interdigital membranes present. Tail about the 
same length as head and body, bicolored, and covered by short, spiny, and clearly visible hairs; 
white hairs present on ventral caudal surface from base to midlength; tail tip with very short 
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FIG. 5. Collection localities of specimens of the Dubosti Group of Neacomys examined in the present study. 
Numbers are keyed to entries in our gazetteer (appendix 1). 


(1 mm) terminal tuft; caudal scales small, arranged in annular series; each caudal scale with 
three subequal hairs inserted along its posterior margin. 

Skull small and delicate in dorsal view (fig. 7), with straight anterior nasal margins, notably 
broad rostrum, and shallow zygomatic notches; posterior nasal terminus usually slightly 
pointed, extending beyond maxillary-frontal suture; premaxillaries terminating anterior to 
nasals; lacrimal bone small and visible in dorsal view, usually in broad contact with frontal 
bones; supraorbital margins slightly convergent anteriorly; interorbital region relatively broad; 
supraorbital beads developed as projecting shelves; lateral expansion of the parietal restricted 
to the dorsal cranial surface. Incisive foramina small and usually teardrop shaped, not extend- 
ing posteriorly to level of M1s; maxillary portion of incisive septum (dividing the left and right 
foramina) narrow. Zygomatic plate narrow. Palate with two posterolateral pits on each side. 
Auditory bullae small and globular, with short and narrow eustachian tubes; periotic bone 
extends anteriorly to internal carotid canal (except MPEG 40435, MPEG 40443, and MPEG 
40439, in which the periotic does not reach the internal carotid canal), but does not enter the 
canal. Subsquamosal fenestra large (almost half the size of the postglenoid foramen); hamular 
process of the squamosal long. Paraoccipital process narrow and small and separated from the 
auditory bullae. Sphenopalatine foramen large; alisphenoid strut absent; carotid circulation 
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FIG. 6. Habitat of Neacomys marajoara at Tauari Farm, municipality of Chaves, Marajé Island, state of Para, 
Brazil, Photo by Rogério V. Rossi, 


pattern usually derived’ (pattern 2, as identified by the presence of a large stapedial foramen 
and a large posterior opening of the alisphenoid canal, and by the absence of a squamosal- 
alisphenoid groove and a sphenofrontal foramen; Voss, 1988). 

First upper molar (M1) usually with broad, flat, and undivided anterocone; anteroloph 
usually fused with anterolabial conule (such that the anteroflexus is not distinguishable); M3 
small; labial cusps (paracone and metacone) usually smaller than lingual cusps (protocone and 
hypocone); m1 anterocone undivided. 

Mandible with mental foramen opening laterally; capsular process of lower incisor alveolus 
present, but indistinct (reduced to a slight rounded elevation), approximately at same height 
as coronoid process. 

KARYOTYPES: Karyotypes of specimens from the type locality were obtained by Silva et al. 
(2015), and karyotypes from specimens collected elsewhere on Marajé Island were obtained 


” Most specimens of Neacomys marajoara (seven out of nine) exhibit carotid circulation pattern 2, but one 
specimen (MPEG 40446) exhibits circulatory pattern 1, and another (MPEG 40431) exhibits both patterns 
(on opposite sides of the skull). 
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FIG. 7. Dorsal, ventral, and lateral cranial views and lateral view of mandible of Neacomys marajoara (MPEG 
40432, holotype). 


by Oliveira da Silva et al. (2019), who described a chromosomal complement of 2” = 58 anda 
fundamental number (FN) = 70 (table 5). 

TAXONOMIC COMPARISONS: Neacomys marajoara differs from N. dubosti in dorsal pelage color 
(dark brown finely sprinkled with orange versus light to dark brown finely sprinkled with orange 
in N. dubosti), and by its bicolored tail (the tail is usually unicolored in N. dubosti), narrower inter- 
orbital region, teardrop-shaped incisive foramina (the lateral margins of the incisive foramina are 
usually subparallel in N. dubosti), globular auditory bullae (the bullae are usually flask-shaped in N. 
dubosti), and carotid circulation usually pattern 2 (versus pattern 1 in N. dubosti). 


2020 
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TABLE 4. Geographic and morphological comparisons among members of the Dubosti Group of Neacomys. 


Distribution 


Dorsal 
pelage 


Ventral 
pelage 


Tail 


Nasal bones 


Interorbital region 


Paraoccipital process 


Subsquamosal fenes- 
trae 


Sphenopalatine fora- 
men 


Lateral expansion of 
parietal 


Carotid circulation 
Incisive foramina 


Maxillary septum of 
incisive foramina 


Zygomatic plate 


Auditory bullae 


M1 anteroloph 


dubosti 


NE Amazonia 


light to dark brown 
finely sprinkled with 
orange 


pure white to buffy 
white 


usually unicolored 


broad 


restricted 


pattern 1 


small and usually 
subparallel 


narrow 


usually flask-shaped 


xingu 


SE Amazonia 


orange-brown 
sprinkled with black 


pure white to buffy 
white 


usually unicolored 
expanded anteriorly 


narrow 


usually separated 
from auditory bullae 


small 


small 


slightly expanded to 
restricted 


usually pattern 14 
small and subparallel 
wide 

narrow to broad 


usually globular 


fused 


marajoara 
SE Amazonia 


dark brown 
finely sprinkled with 
orange 


pure white to yellow- 
ish white 


bicolored 
straight anteriorly 


relatively broad 


separated from audi- 
tory bullae 


large 


large 


restricted 


usually pattern 2° 


small and usually 
teardrop shaped 


narrow 


narrow 


globular 


fused or distinct 


* We found variation in this character, please see Neacomys xingu description. 
> We found variation in this character, please see Neacomys marajoara description. 


vossi 


SE Amazonia 


brown 
sprinkled with 
orange and/or black 


pure white to buffy 
white 


usually slightly bicol- 
ored 


expanded anteriorly 


narrow 


close to the auditory 
bullae 


usually small 
large 


slightly expanded to 
restricted 


pattern | 


small and subparallel 


usually wide 


narrow 


globular 


usually distinct 
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Neacomys marajoara differs from N. xingu in dorsal pelage color (dark brown finely 
sprinkled with orange versus orange-brown sprinkled with black in N. xingu), and by its 
bicolored tail (the tail is usually unicolored in N. xingu), anteriorly straight (versus anteriorly 
expanded) nasal bones, broader interorbital region, narrower maxillary septum of the inci- 
sive foramina, larger subsquamosal fenestrae, and carotid circulation usually pattern 2 (ver- 
sus usually pattern 1). 


26 AMERICAN MUSEUM NOVITATES NO. 3958 


Neacomys marajoara differs from N. vossi in dorsal pelage color (dark brown finely sprin- 
kled with orange versus brown sprinkled with orange and/or black in N. vossi), and by its 
anteriorly straight (versus anteriorly expanded) nasal bones, broader interorbital region, para- 
occipital processes separated from the auditory bullae (versus processes close to the auditory 
bullae), larger subsquamosal fenestrae, carotid circulation usually pattern 2 (versus pattern 1), 
and narrower maxillary septum of the incisive foramina. 

Karyotypically, Neacomys marajoara differs from other species in the Dubosti Group by 
having a diploid chromosomal complement of 58 (versus 2n = 64 in N. dubosti; table 5), a 
uniquely large FN of 70 autosomal arms (versus FN = 64-68 in other group members), and 
submetacentric sex chromosomes (at least one of the sex chromosomes is acrocentric in N, 
dubosti and “species 2”), 

DisTRIBUTION: Neacomys marajoara has been recorded only on Marajé Island, state of 
Para, Brazil (fig. 5). 

Erymo.ocy: The specific epithet marajoara is to be treated as a noun in apposition, It is 
derived from Tupi Guarani (a language family that comprises many different indigenous Ama- 
zonian dialects) and denotes a native of Marajé Island (the type locality). 

FreLp Notes: Of the nine specimens for which trapping information is available, eight 
were captured in pitfall traps and one was taken in a Sherman trap placed on the ground. One 
paratype (MPEG 40435) was pregnant with three fetuses. 

REMARKS: Neacomys marajoara was previously reported in the literature as “Neacomys sp.’ 
(in part) by Silva at al. (2015) and as “Neacomys sp. D” by Oliveira da Silva et al. (2019). 


Neacomys vossi, new species 


Vosss Spiny Mouse 
Figures 8, 9 


Ho.otyPe: The holotype (UFPA 1583) is an adult female (age class 2) collected on 2 April 
2013, by Ana Cristina Mendes Oliveira (original field number JB 031). The specimen consists 
of a stuffed skin, skull, and skeleton, all in good condition; additionally, a tissue sample is pre- 
served in ethanol, and a partial cytochrome b sequence that we obtained from it has been 
deposited in Genbank with accession number MT462028. 

Type Loca.irty: Boca do Rato, on the right bank of the Rio Tapajos, Itaituba municipality, 
state of Para, Brazil (5°14’S, 56°56’W, fig. 5). 

Dracnosis: Neacomys vossi is a small species (table 3) that differs from congeneric taxa by 
the following combination of craniodental traits: skull delicate; interorbital region narrow; 
nasal bones expanded anteriorly; supraorbital margins convergent anteriorly; subsquamosal 
fenestra usually small (about %4 the size of the postglenoid foramen on each side of the skull); 
paroccipital processes close to the auditory bullae; sphenopalatine foramen large; carotid cir- 
culation pattern 1 (sensu Voss, 1988); maxillary part of incisive septum (between the incisive 
foramina) usually wide; M1 usually with slightly flat, narrow, and undivided anterocone; and 
M1 anteroloph usually distinct (not fused with the anterolabial conule). 
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MoRPHOLOGICAL DEscriPTiION: Dorsal pelage brown sprinkled with orange and/or black 
(fig. 8); ventral pelage varying from pure white to buffy white; very thin orange lateral line. Super- 
ciliary, genal, and mystacial vibrissae blackish and long (extending behind ears when laid back 
alongside the head); submental vibrissae absent; interramal vibrissae short and white, Ears small 
and rounded; post-auricular hairs gray-based with orange tips, forming an orange tuft behind 
each pinna. Ungual tufts white, longer than claws; fore- and hind feet covered dorsally with buffy- 
cream hairs; hind feet narrow and elongate with small interdigital membranes. Tail about the 
same length as head and body, usually slightly bicolored (except in UFPA 1277, 1417, 1444, 1654, 
and 1736 which have unicolored tail) and covered by spiny and clearly visible hairs; tail tip of the 
tail with very short (approximately 1.4 mm long) terminal tuft; caudal scales small, arranged in 
annular series; each caudal scale with three subequal hairs inserted along its posterior margin. 

Skull small and delicate in dorsal view (fig. 9); with expanded anterior nasal margins; 
notably broad rostrum and shallow zygomatic notches; posterior nasal terminus usually 
flat, extending beyond the maxillary-frontal suture; premaxillaries terminating slightly 
anterior to nasals; lacrimal bone small and visible in dorsal view, equally contacting maxil- 
lary and frontal bones; supraorbital margins convergent anteriorly; interorbital region nar- 
row; supraorbital beads developed as projecting shelves; lateral expansion of the parietal 
restricted to the dorsal surface (except in UFPA 1227 and UFPA 1577, in which the parietal 
is slightly expanded ventrally near the squamosal root of the zygomatic arch), Incisive 
foramina small with subparallel lateral margins, not extending posteriorly to the level of 
M1s; maxillary portion of incisive septum (dividing the left and right foramina) usually 
wide. Zygomatic plate narrow. Palate with two posterolateral pits on each side of the palate. 
Auditory bullae small and globular, with short and narrow eustachian tubes; periotic bone 
extending anteriorly to internal carotid canal, but usually not entering into it (except in 
UFPA 1227, 1417, 1530, 1654, 1417, and 1736, in which the periotic does enter the internal 
carotid canal). Subsquamosal fenestra usually small (about % the size of the postglenoid 
foramen); hamular process of the squamosal long. Paraoccipital process narrow, small, and 
close to the auditory bullae (Sanchez-Vendizu et al., 2018: fig. 3C). Sphenopalatine fora- 
men large; alisphenoid strut absent; carotid circulation primitive (pattern 1, as identified 
by retaining a well-developed squamosal-alisphenoid groove and sphenofrontal foramen, 
both indicative of the presence of the supraorbital branch of the stapedial artery; Voss et 
al., 1988). 

First upper molar (M1) usually with slightly flat and undivided anterocone; anteroloph 
usually distinct (not fused with the anterolabial conule); M3 small; labial cusps (paracone, 
metacone) usually taller than lingual cusps (protocone, hypocone); ml anteroconid 
undivided. 

Mandible with mental foramen opening laterally; capsular process of lower incisor alveolus 
present, but indistinct (reduced as a slight, rounded elevation), approximately at same height 
as coronoid process. 

TAXONOMIC COMPARISONS: Neacomys vossi differs from N. dubosti in dorsal pelage color 
(brown sprinkled with orange and/or black versus dark brown finely sprinkled with orange in 
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FIG. 8. Dorsal and ventral views of the skin of Neacomys vossi (UFPA 1583, holotype). 
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FIG. 9. Dorsal, ventral, lateral cranial views and lateral view of mandible of Neacomys vossi (UFPA 1583, 
holotype). 
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_ N. dubosti), and by its narrower interorbital 
region, and globular auditory bullae (the bullae 
are usually flask shaped in N. dubosti). 

Neacomys vossi differs from N. xingu in dor- 
sal pelage color (brown sprinkled with orange 
and/or black versus orange-brown sprinkled 
with black in N. xingu), larger sphenopalatine 
foramen, and an M1 anteroloph that is usually 
distinct (versus fused in N. xingu). 

Karyotypically, Neacomys vossi differs from 
other species in the Dubosti Group by having a 
Bee chromosome complement of 58 (versus 
2n = 64 in N. dubosti and 2n = 64 in “species 2”), 
an FN of 68 autosomal arms (versus FN = 70 in 
N. marajoara, FN = 64 in N. xingu, and FN = 66 
in “species 2”), and submetacentric sex chromo- 
somes (at least one of the sex chromosomes is 
acrocentric in N. dubosti and “species 2”). 

DISTRIBUTION: Neacomys vossi has been col- 
lected on the right bank of the upper and middle 
Tapajos River and on the left bank of lower 
Xingu (fig. 5). According to our records, the spe- 
FIG. 10. Robert S. Voss in the field (Venezuela, ‘ies appears to be endemic to the Tapajos center 
1979; photo by Paul Kaarakka). of endemism (Silva et al., 2005). 

ETyMoLoGy: Named in honor of Robert S. 
Voss (fig. 10), curator of mammals at the American Museum of Natural History, New York, for 
his extensive contributions to our knowledge of Neotropical mammals, especially the taxonomy 
of Neacomys in northeastern Amazonia. 

FiELD Notes: Among the specimens we examined for which trapping information is avail- 
able, all 17 were captured in pitfall traps, two (UFPA 1277 and 1284) in upland (terra firme) 
forests and three (UFPA 1444, 1487, and 1583) in primary forest of unknown character. 

REMARKS: Neacomys vossi was previously reported in the literature as “Neacomys sp. A’ by 
Oliveira da Silva et al. (2017, 2019). 





Neacomys xingu, new species 


Xingu Spiny Mouse 
Figures 11, 12 


Ho.otyPe: The holotype (UFMT 1268) is an adult female (age class 5), collected on 28 
August 2009 by Cleuton Lima Miranda (original field number PSA 242) in a pitfall trap. The 
specimen is preserved as a skin, skull, and skeleton in good condition; additionally, a tissue is 
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preserved in alcohol, and a partial cytochrome b sequence that we obtained from it has been 
deposited in Genbank with accession number MT462060. 

Type Loca.ity: Flona Tapirapé-Aquiri, Maraba, state of Para, Brazil (5°46'S, 50°32'W, 
fig. 5). 

Diacnosis: Neacomys xingu is a small species (table 3) that differs from congeneric taxa 
by the following combination of craniodental traits: skull delicate; interorbital region narrow; 
nasal bones expanded anteriorly; supraorbital margins convergent anteriorly; subsquamosal 
fenestra small (about % the size of the postglenoid foramen on each side of the skull); paraoc- 
cipital process separated from the auditory bullae; sphenopalatine foramen small; carotid cir- 
culation usually pattern 1 (sensu Voss, 1988); maxillary part of incisive septum (between the 
incisive foramina) wide; M1 usually with flat and undivided anterocone; and M1 anteroloph 
fused with the anterolabial conule. 

MorRPHOLOGICAL DescriPTION: Dorsal pelage orange-brown, sprinkled with black (fig. 
11); ventral pelage varying from pure white to buffy white and separated from the dorsal pelage 
by a thin orange lateral line, Superciliary, genal, and mystacial vibrissae blackish and long 
(extending behind ears when laid back alongside the head); submental vibrissae absent; inter- 
ramal vibrissae short and white; Ears small and rounded; postauricular hairs gray based with 
orange tips, forming an orange tuft behind each pinna, Ungual tufts white, longer than claws 
in length; fore- and hind feet covered dorsally with buffy-cream hairs; hind feet narrow and 
elongate with small interdigital membranes present. Tail about the same length as head and 
body, usually unicolored (except in PSA 069, MPEG 39901, and MPEG 42019, in which the 
tail is dark above and paler below), and covered by short, spiny, and clearly visible hairs; tail 
tip with very short (1-2 mm) terminal tuft; caudal scales small, arranged in annular series; each 
caudal scale with three subequal hairs inserted along its posterior margin, 

Skull small and delicate in dorsal view (fig. 12); with anteriorly expanded nasal margins; 
notably broad rostrum and shallow zygomatic notches; posterior nasal terminus slightly 
pointed, extending beyond the maxillary-frontal suture; premaxillaries terminating slightly 
anterior to nasals; lacrimal bone small and visible in dorsal view, equally contacting the maxil- 
lary and frontal bones; supraorbital margins convergent anteriorly; interorbital region narrow; 
supraorbital beads developed as projecting shelves; lateral expansion of the parietal restricted 
to the dorsal cranial surface (except in MCN-M 1404 and MPEG 39901, in which the parietal 
is slightly expanded ventrally near the squamosal root of the zygomatic arch). Incisive foramina 
small with subparallel lateral margins, not extending posteriorly to level of M1s; maxillary 
portion of incisive septum (dividing the left and right foramina) usually wide. Zygomatic plate 
varying from broad to narrow. Palate with two posterolateral pits on each side. Auditory bullae 
small and usually globular, with short and narrow eustachian tubes; periotic bone extends 
anteriorly to the internal carotid canal but does not enter it (except in MPEG 42715 and 
MCN-M 1404, in which the periotic does not reach the internal carotid canal). Subsquamosal 
fenestra small (about % the size of the postglenoid foramen); hamular process of the squamosal 
long. Paraoccipital process narrow and small and usually separated from the auditory bullae 
(Sanchez-Vendizu et al., 2018: fig. 3C). Sphenopalatine foramen small (except in MPEG 41991, 
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10 mm 


FIG. 11. Dorsal and ventral views of the skin of Neacomys xingu (UFMT 1268, holotype). 
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whose sphenopalatine foramen is large); alisphenoid strut absent; carotid circulation pattern 
usually primitive’ (pattern 1, as identified by retaining a well-developed squamosal-alisphenoid 
groove and sphenofrontal foramen, both indicative of the presence of the supraorbital branch 
of the stapedial artery; Voss, 1988: fig. 18A, B), 

First upper molar (M1) usually with slightly flat and undivided anterocone; anteroloph 
usually fused with the anterolabial conule (such that the anteroflexus is not distinguishable); 
M3 small; labial cusps (paracone, metacone) usually taller than lingual cusps (protocone, hypo- 
cone); ml anteroconid undivided. 

Mandible with mental foramen opening laterally; capsular process of lower incisor alveolus 
present, but indistinct (reduced to a slight rounded elevation), approximately at same height 
of coronoid process, 

TAXONOMIC Comparisons: Neacomys xingu differs from N. dubosti in dorsal pelage color 
(orange-brown sprinkled with black versus light to dark brown finely sprinkled with orange), 
and by its narrower interorbital region, globular auditory bullae (the bullae are usually flask 
shaped in N, dubosti), and carotid circulation usually pattern 1 (versus always pattern 1 in N, 
dubosti). 

Karyotypically, Neacomys xingu differs from all other members of the Dubosti Group by 
having a diploid chromosomal complement of 58 (versus 2n = 64 in N. dubosti and 2n = 54 in 
“species 2”), a uniquely small FN of 64 autosomal arms (versus FN = 66-60 in other species), 
and submetacentric sex chromosomes (at least one sex chromosome is acrocentric in N. dubosti 
and “species 2”). 

DIsTRIBUTION: Neacomys xingu has been collected on the right bank of the lower Xingu 
River and in the region of Serra de Carajas, in southeastern Para state (fig. 5), According to 
cytogenetic data (Di-Nizo et al., 2017; Oliveira da Silva et al., 2019), the species also ranges 
southward into Vila Rica in northeastern Mato Grosso state. These records suggest that the 
species is restricted to the Xingu center of endemism (Silva et al., 2005). 

Erymo.ocy: The specific epithet xingu is to be treated as a noun in apposition. The species 
name refers to the Xingu center of endemism, delimited by the Xingu and Tapajds rivers, where 
the species occurs. 

REMARKS: Neacomys xingu was previously reported in the literature as “Neacomys clade 7” 
by Patton et al. (2000), “Neacomys sp.” (in part) by Silva at al. (2015), “Neacomys sp.” by Di-Nizo 
et al. (2017) and Brandao et al, (2019), and “Neacomys sp. C” by Oliveira da Silva et al. (2019). 


DISCUSSION 


This study includes the largest mitochondrial DNA dataset yet assembled for Neacomys, 
the results of which support generic monophyly, as have all previous analyses of Cyth 
sequence data (Catzflies and Tilak, 2009; Silva et al., 2015; Hurtado and Pacheco, 2017; 
Sanchez-Vendizu et al., 2018). Our analyses additionally reveal the existence of 17 distinct 


* Most specimens of Neacomys xingu haye carotid circulation pattern 1, whereas four specimens had circula- 
tory pattern 2 (MPEG 40446, MPEG 39901, MCN-M 1403, and MCN-M 1404), 
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FIG. 12. Dorsal, ventral, lateral cranial views and lateral view of mandible of Neacomys xingu (UFMT 1268, 
holotype). 


lineages clustered in four main groups. Two of these groups correspond to the Paracou and 
Spinosus groups of Hurtado and Pacheco (2017) and Sanchez-Vendizu et al. (2018), but the 
inclusion of new sequence data from southeastern Amazonia in our analyses make it expedi- 
ent to restrict the membership of their Tenuipes Group and to recognize a fourth clade, 
which we refer to as the Dubosti Group. 

The three new species herein described—N. marajoara, N. xingu, and N. vossi—belong to 
the Dubosti Group. These species, together with an unnamed member species (“species 2”), 
have genetic divergences varying from 3.9% to 9.4%, comparable to distances previously 
reported for recently described species such as N. vargasllosai and N. spinosus (7.8%; Hurtado 
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TABLE 5. Karyotypes * of Neacomys species in the Dubosti Group. 








2n FN gA gB gC sD Xx Y Source 
dubosti 64 68 — 6 — 56 SM A Silva et al., 2015 
marajoara 58 70 _ 14/8 = — 42/48 SM SM __ Silva et al., 2015 
vossi 58 68 -- 12 ~ 44 SM SM Oliveira da Silva et al., 2017 
xingu 58 64 - 8 _ 48 SM SM _ Silva et al., 2015 
species 2 54 66 8 6 ~ 38 A A Oliveira da Silva et al., 2017 


4 Diploid number (2); fundamental number (FN); group A (gA) = large metacentric and submetacentric autosomes; 
group B (gB) = medium and small metacentric or submetacentric autosomes; group C (gC) = medium and small sub- 
telocentric autosomes; group D (gD) = large, medium, and small acrocentric autosomes, Abbreviations for sex chromo- 
some (X and Y) morphology: M = metacentric, SM = submetacentric, ST = subtelocentric, A = acrocentric. 


and Pacheco, 2017), and N. macedoruizi and the “upriver” clade of N. minutus (4.9%; Sanchez- 
Vendizti et al., 2018). Although N. marajoara and N. xingu have the lowest genetic divergence 
value (3.9%) reported for any pair of Neacomys species, this value is within the known range 
of interspecific mammalian distances (Baker and Bradley (2006). Additionally, cytogenetic 
studies have revealed karyotypic differences among these lineages (Silva et al., 2015; Oliveira 
da Silva et al., 2017, 2019). These, along with qualitative-morphological data reported in our 
study, provide evidence of reproductive isolation and nuclear-gene divergence, respectively. 

Carotid circulation patterns have long been considered a reliable character for diagnosing 
sigmodontine species (Voss, 1988; Carleton and Musser, 1989; Voss and Carleton, 1993; Step- 
pan, 1995; Carleton and Olson, 1999; Weksler, 2006; Hurtado and Pacheco, 2017; and Sanchez- 
Vendizu et al., 2018). However, intraspecific polymorphism in this character was reported by 
Voss (1991), and our observations provide additional evidence that oryzomyine species can be 
polymorphic for carotid circulatory traits. Carotid polymorphisms in N. marajoara and N. 
xingu may be related to the recency of speciation events from which these lineages arose, as 
evidenced by their minimal genetic divergence, Because these are not the only congeneric spe- 
cies that exhibit carotid polymorphisms (T.B.ES,, unpublished), observations of sample differ- 
ences in carotid traits should be interpreted cautiously unless accompanied by other evidence 
of taxonomic divergence. 

The discovery of new species of Neacomys in southeastern Amazonia is not unexpected, 
in view of similar recent discoveries in western and northeastern Amazonia (Patton et al., 
2000; Voss et al., 2001; Hurtado and Pacheco, 2017; and Sanchez-Vendizu et al., 2018) and 
previous reports of karyotypic variation and genetic divergence among southeastern forms 
(Patton et al., 2000; Catzeflis and Tilak, 2009; Silva et al., 2015; Oliveira da Silva et al., 2017; 
Hurtado and Pacheco, 2017; and Sanchez-Vendizt et al., 2018). Although southeastern Ama- 
zonia has long been accessible to collectors, our study underscores the fact that faunal inven- 
tories and taxonomic studies are still needed in this region, which is among the most 
extensively deforested and environmentally threatened parts of the Amazonian biome, due 
to industrial agriculture, road construction, mining, and hydroelectric dams (Fearnside, 
1989, 2005; Gascon et al., 2001). 
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APPENDIX 1 


SPECIMENS EXAMINED 

The localities of specimens from the Dubosti Group that we examined and from which we 
analyzed sequence data are listed below. The largest administrative units within each country 
(state, department, or province) are italicized. Verbatim locality names are placed in quotes, 
but others have been emended for brevity and to conform with current geographic usage. For 
each locality we provide geographic coordinates and the collection number of the voucher 
specimens examined. Numbers identify locality symbols plotted on map (fig. 5). The paratypes 
of our new species are designated here. 


BRAZIL 

1, Amapa, Serra do Navio km 190 efa. (00°59'S, 52°03’ W; USNM 461589; N. dubosti), 

2. Amapd, Rio Amapari (00°59’S, 52°03’°W; USNM 461573, 461587; N. dubosti). 

3. Pard, Almeirim Municipality, Reserva Biologica Maicuru (00°50°N, 53°56’W; CN 240; N. 
dubosti). 

4, Para, Chaves municipality, Marajé Island, Tauari Farm (00°39’S, 50°11°W; MPEG 40429, 
40431, 40434, 40435, 40439-40441, 40443, 40446 [paratypes]; MPEG 40432 [holotype]; N. 
marajoara). 

5. Para, Mangabal community on right bank of Rio Tapajds, Itaituba (5°35’S, 57°7°W; UFPA 
1736 [paratype]; N. vossi). 

6. Para, Itaituba Municipality, Boca do Rato on right bank of Rio Tapajés (6°7’S, 57°36 W; 
UFPA 1391, 1520, 1577 [paratypes]; UFPA 1583 [holotype]; N. vossi). 

7. Para, Ttaituba, Mamae-ana (Canta-galo) on right bank of Rio Tapajos (5°48’S, 57°24’ W; 
UFPA 1277, 1284, 1487, 1647, 1654 [paratypes]; N. vossi). 

8, Parad, Itaituba (5°40’S, 57°00’ W; UFPA 1467, JB 27; N. vossi). 

9. Pard, Itaituba, Sao Martins on right bank of Rio Tapajos (6°07’S, 57°36’W; UFPA 1691 
[paratype]; N. vossi). 

10, Pard, Itaituba, Jacaré on right bank of Rio Tapajés (5°41’S, 57°14" W; UFPA 1444 [paratype]; 
N. vossi). 

11. Para, 18 km south and 19 km west of Altamira- Agrovila da Uniao (03°22°S, 52°23'W; 
USNM 521468, 521539; N. vossi), 

12. Para, Itaituba, Br 230, Itaituba-Jacareacanga Km 212 on Transamazonica Highway (5°40'S, 
56°45°W; USNM 545957-545959, 545960-545964; N. vossi). 

13. Para, Maraba Municipality, Flona Tapirapé-Aquiri (5°51°S, 50°31’W; UFMT 1273, 1275 
[paratypes]; MPEG 39901, 41804, 41805, 41991, 41996, 42019; MPEG-PSA 46, 69 [para- 
types]; UFMT 1268: [holotype]; N. xirgu). 

14. Pard, East bank Rio Xingu 52 km SSW Altamira (3°39°S, 52°227W; USNM 549553; N. 
xingu). 

15. Pard, 54 km south and 150 km west of Altamira (3°41’S, 51°45°W; USNM 519180, 521464- 
521467; N. xingu). 
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16. Pard, Flona Carajas, Parauapebas (6°3’S, 50°35°W; MCN-M 1401-1404, 1443 [paratypes]; 
N. xingu). 

17. Para, left bank of Rio Tapajés (5°17’S, 57°8°W; MPEG 42838, 42901; Neacomys sp. 2). 

18. Para, Juruti (2°09°S, 56°5°W; JUR 19, 42, 82; Neacomys sp. 2). 

19. Para, Penedo on left bank of Rio Tapajdés (5°30°S, 57°6' W; UFPA 1227, 1413, 1417, 1530; 
Neacomys sp. 2). 

20. Pard, Reserva Bioldgica Maicuru. 00°50°N. 53°56’W (CN 240: Neacomys dubosti). 


FRENCH GUIANA 

21. Paracou near Sinnamary (05°17°N. 52°55’W; AMNH 267569 [holotype]; N. dubosti) 
22. Cacao, Roura, Cayenne (04°34’N, 52°27°W; MHNG 1963.025: N. dubosti). 

23. Route de Kaw, Cayenne (04°36’N, 52°15’°W; MHNG 1894.013: N. dubosti). 


SURINAM 

24. Nickerie District, Sipaliwini Airstrip (02°02’N, 56°07'W; CM 76846: N. dubosti). 

25. Sipaliwini, Werehpai Camp (02°23’N, 56°42’W: ROM 120726: N. dubosti). 

26. Sipaliwini, Iconja Landing on Sipaliwini River (01°59’N, 56°05’ W: ROM 120288: N. dubosti). 

27. Marowijne, Oelemarie (03°06'N, 54°31’ W; CM 76835-76837, 76839, 76840, 76842-76844; 
N. dubosti). 
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APPENDIX 2 


PHYLOGENETIC RELATIONSHIPS OF NEACOMYS BASED ON MAXIMUM-LIKELIHOOD ANALYSIS 


(For simplicity, clades shared with our Bayesian results [in fig. 2A—C] are 
cartooned. Branch support values are bootstrap frequencies.) 
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APPENDIX 3 


RESULTS OF PRINCIPAL COMPONENTS ANALYSES OF SPECIES 
IN THE DUBOSTI GROUP OF NEACOMYS 
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